Interfacial phase transitions in twinning-plane superconductors 
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Abstract 

Within the framework of Ginzburg-Landau theory we study the rich variety of interfacial phase 

G 

transitions in twinning-plane superconductors. We show that the phase behaviour strongly de- 
■ pends on the transparency of the twinning plane for electrons measured by means of the coupling 

parameter «tp • By analyzing the solutions of the Ginzburg-Landau equations in the limit of per- 
fectly transparent twinning planes, we predict a first-order interface derealization transition for 
all type-I materials. We further perform a detailed study of the other limit in which the twinning 
plane is opaque. The phase diagram proves to be very rich and fundamentally different from the 

o 

transparent case, recovering many of the results for a system with an external surface. In particular 
both first-order and critical derealization transitions are found to be possible, accompanied by a 
■^j- . first-order depinning transition. We provide a comparison with experimental results and discuss 

^NJ , the relevance of our findings for type-II materials. 

o 

PACS numbers: 05.70.Np , 74.25.Dw , 61.72.Mm , 64.60.Fr 

O . 
-t— > ' 

s 

£ 

O 

o 

> 
X 



*Present address: Department of Mathematics, Heriot-Watt University, Edinburgh EH14 4AS, United King- 
dom. 



1 



I. INTRODUCTION 



In recent years, local enhancement of superconductivity has been predicted to provide 
the mechanism to induce several intriguing interfacial phase transitions in type-I supercon- 
ductors [Q, 0, ||, f|, |[. Typical phase diagrams are calculated using the Ginzburg-Landau 
(GL) theory in which the enhancement is accounted for by allowing the extrapolation length 
b to be negative. The microscopic origin of this parameter remains an unsolved problem 
making the experimental verification of the theoretical results non-trivial. So far, the most 
feasible realization of a negative extrapolation length seems to originate from the concept of 
twinning-plane superconductivity (TPS), a well understood phenomenon that occurs, e.g., 
in Sn, In, Nb, Re and Tl 0. A twinning plane (TP) is a defect plane representing the bound- 
ary between two single-crystal regions or twins and, consequently, the physics encountered 
in the behaviour of a superconducting/normal interface near TP's is the natural analogy 
to grain-boundary wetting or interface depinning, a topic which has been well studied in 
magnetic systems [0, f| 0. 

The characteristic feature of the original GL approach of TPS is the a priori assumption 
that the TP is perfectly transparent for electrons at the microscopic level which implies 
that the superconducting order parameter if) is continuous at the TP ||. Subsequent ex- 
tensions of the theory relax this assumption allowing a discontinuity in ip [K], [11], [1^, |13|| . 
More specifically, a second phenomenological parameter, a TP , is introduced to describe the 
coupling between the twins such that, by means of cktp, one can mimic the effect of micro- 
scopically tuning the TP from completely transparent to completely opaque for electrons. 
In this paper we present an overview of the variety of interfacial phase transitions in the two 
limiting cases to develop a thorough understanding of the influence of the transparency. 

Earlier work j3J has focused on the case of mixed bulk boundary conditions with the 
bulk normal (N) phase on one side and the bulk superconducting (SC) phase on the other 
side of the TP, at bulk two-phase coexistence. This is appropriate for the study of the 
proper depinning transition of an interface that is initially pinned at the TP. Here we choose 
to settle for the configuration of equal bulk conditions, that is, we impose the bulk N 
phase on both sides of the TP. In so doing we are no longer restricted to the case of bulk 
two-phase coexistence and this allows us to establish the complete magnetic field versus 
temperature phase diagram for a given material. This type of diagram is accessible to 
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experimental verification and is relevant for comparing the present results with known TPS 



phase diagrams p, 12, 13 



The solutions of the GL equations strongly depend on the boundary conditions imposed 
at the TP itself which in turn relates to the level of transparency, i.e. the value of otp- 
For highly transparent planes, corresponding to the limit «tp — > 0, it is natural to consider 
fully symmetric profiles for the order parameter. In the opposite limit, «tp — ► oo, the 
TP is completely opaque for electrons and both sides are largely independent. In this case 
there is a wide range of possible solutions, including profiles with ip identically zero at one 
side of the TP. The latter are refered to as wall solutions, since they are equivalent to the 
ones found in a type-I superconductor with an external surface or wall characterized by a 
negative extrapolation length b ||. Therefore we anticipate that, in the opaque limit, 
we will recover to a great extent the results of a wall system. This is very different from 
the case of complete transparency, for which drastic qualitative modifications are predicted 
compared to the case with a wall. 

The outline of the paper is as follows. In the next section we collect the main ideas of 
the GL theory applied to twinning-plane superconductors. Section [Til] covers the results 
for perfectly transparent TP's. We calculate in detail the phase diagrams and provide a 
comparison with the predicted TPS diagrams as described in Ref. ||. The fully opaque 
system is the subject of Section ffVl. We present a classification of the various solutions and 
establish their stability to derive the phase behaviour. We summarize our main results and 
discuss the experimental relevance in Section [V|. 



II. GINZBURG-LANDAU THEORY FOR TWINNING-PLANE SUPERCON- 
DUCTORS 

We consider a type-I superconductor with a TP located at x = and impose on both 
sides the N phase, with ip = 0, as the bulk condition. The GL free-energy functional has 
the form 

/+oo 
£#,A]dx + r T p(?/>-, (1) 
-oo 

with the free-energy density Q given by 



As usual, e oc T — T c , where T c is the bulk critical temperature which must be distinguished 
from the second critical temperature in the system, T Cj tp, below which local superconduc- 
tivity sets in at the TP in zero magnetic field. Since T c tp was experimentally || proved to 
be only slightly higher than T c , the use of the GL theory is justified. Further, (3 > is a 
stabilizing parameter and A is the vector potential. We choose the applied magnetic field 
H = He z parallel to the TP. Using the notation ?/>_ = an d ?p+ = ip{0 + ), the local 

contribution Ttp in ([]]) reads 

rxpW-.W = ^W+l 2 + + ^ |* + - (3) 

The first term, with b < 0, describes the enhancement of superconductivity and was intro- 
duced by Khlyustikov and Buzdin || to reproduce theoretically the observed TPS phase dia- 
grams. The phenomenological parameter b is the extrapolation length and can be related to 
the temperature difference T c — T Ci tp- In addition, we have followed others @, |10 , TT|, |12, [13 



by adding a second term in @ to describe the coupling between the twins. In so doing, we 
allow the SC wave function to be discontinuous across the TP, hence in general ^_ ^ 
The coupling constant q;tp can be expressed in terms of the Fermi velocity and either the 
transmission or reflection coefficient for electrons, thus fully in terms of microscopic prop- 



erties ||11|| . We note that for «tp > 0, the phase of the wave function is continuous at the 
TP, while for axp < a phase jump of n can occur |I0| . We omit the latter possibility and 
restrict our attention to «tp > 0. 

In what follows we assume translational invariance in the y- and z-directions and choose 
the gauge so that A = (0, A(x),0). It proves to be convenient to adopt the rescaling 
introduced in earlier work |2j using the two basic length scales of the superconductor, i.e. 
the zero-field coherence length £ and the magnetic penetration depth A defined by 



^2 = ^ A 2 = m P / 4 N 

2m\e\ ' /x ^ 2 1 e 1 



The ratio of A to £ gives the GL parameter k, with n < 1/^/2 for type-I materials. We 
use £ to scale the distances but, for simplicity, retain the notation x for the dimensionless 
coordinate x/£ perpendicular to the TP. For the magnetic quantities A and H we introduce 
the dimensionless a and h defined by 

2eA 2eAVo 
a=—A , h^—f—H, (5) 



and rescale the wave function tp according to 

<p = ip/ipeq, (6) 



where ip eg = J\e\/(3 is the equilibrium value of the SC order parameter for T < T c . Clearly, 
(p attains the value 1 in the bulk SC phase. Finally, we rescale the free energy T divided by 
the surface area S such that 7 = Y(3/ \e 2 ^S), yielding 



l[(p,a] 



J— 00 12k J 

+ + — (<P- -<P + ) 2 - (7) 

b OTP 



The ± refers to the sign of T — T c . Minimization of 7 with respect to (p and a yields the 
well-known GL equations 

(p = ±ip + a 2 ip/n 2 + y? 3 , (8) 

and 

a = acp 2 /K 2 . (9) 

In addition, two coupled boundary conditions are obtained from stationarity with respect 
to and ip + 

ip_ = -t<P- + (<p+-<p-), (10) 

Q?TP 

0+ = + — (V+ - <P-)> (11) 

b axp 

while the vector potential and its first derivative must be continuous at x — 0, 

a(x = 0') = a{x = + ) , a(x = 0") = a(x = + ). (12) 

In the subsequent sections we aim at solving the above equations for both the transparent 
and the opaque limit. 

III. TRANSPARENT TWINNING PLANES 
A. Boundary conditions and solutions 

For highly transparent TP's, «tp - > and we recover the original description of TPS 
with a continuous order parameter at the TP, thus (p + — (p~. Consequently, it is natural to 
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look for fully symmetric solutions for (f(x) of the differential equation (|8|). In practice this 
is done by initially restricting ourselves to one half space, say x > 0. The solution in the 
other half space x < can then be constructed using (f(x) = ip(—x), with <p(x) the solution 
for the semi-infinite system. Note that this makes the problem very similar to the one of a 
semi-infinite system with an external surface |TJ, the only difference emerging from the 
behaviour at x = which is manifested most clearly in the boundary condition for the vector 
potential. Indeed, in the presence of an external surface, stationarity of 7 with respect to 
a(0) results in d(0) = h. In the present situation, however, owing to the assumed symmetry 
in the profile for <f(x) the vector potential a(x) will be anti-symmetric with respect to x, i.e. 
a(— x) = —a(x) and thus obeys the condition a(0) = 0. The boundary conditions for the 
wave function reduce in this limit to 

tp + = -if. = -tp+. (13) 

In the following analysis, we will be interested in two types of solutions distinguished 
by their asymptotic behaviour for \x\ — > 00. The enhancement of superconductivity near 
the TP will typically induce a SC sheath (with the bulk of the system prepared in the N 
phase) and it is precisely the thickness of this surface sheath that characterizes the solution. 
Clearly, if this thickness is finite superconductivity disappears for |x| — ^ 00 and the magnetic 
field penetrates such that 

y?(±oo) = , a(±oo) = h. (14) 
Adapting the terminology of ]l|, H which follows from the analogy to wetting transitions in 



adsorbed fluids |14|, this solution is referred to as a partial wetting state. We remark that 
this class of solutions also includes the so-called null solution without any SC phase in the 
system (thus describing a SC sheath with zero thickness). On the other hand, in the case 
of complete wetting the sheath can be macroscopically thick so that it fully separates the 
N phase from the TP. These solutions are possible only when the magnetic field equals the 
thermodynamic critical field H c at which there is bulk two-phase coexistence between the 
N and SC phase. Instead of Eq.(|14"D they obey the asymptotic conditions 



ip{±oo) = 1 , a(±oo) = 0. (15) 

Remarkably, the calculation of this class of solutions is very simple and analytic results can 
be obtained for the entire type-I regime. Indeed, the asymptotic condition for the vector 



potential along with the boundary condition at x = imply that a = and we are left with 
a single differential equation for (p(x) taking the simplified form, for T < T c , 

(p = -^ + v\ (16) 

with the solution tp(x) = coth[(x + 5)/y/2], for x > 0, and where 5 is determined by the 
boundary condition ([13]) This solution must then be combined with a SC/N interface in 
the limit of \x\ — > oo as necessary to obey the bulk N condition as initially imposed. This 
combination defines a macroscopic SC layer. For general values of k, the SC/N interface, 
as well as the solutions describing a finite SC sheath, must be determined numerically. 
At bulk two-phase coexistence the two solutions described above may occur yielding an 
interface derealization transition (or a wetting transition). This transition describes the 
derealization of the SC/N interface from the TP into the bulk of the material away from 
the TP. We now address the analysis of this phase transition. 

B. Phase diagram at bulk coexistence 

In FigJI] we show the phase diagram for a type-I superconductor with a transparent TP as 
a function of the parameters k and £/&. The magnetic field h is fixed to its coexistence value 
h c = l/v2 |15| . Varying the temperature for a given material, i.e. for fixed k, corresponds 
to traveling horizontally in the diagram. Moving along coexistence towards T c corresponds 
to decreasing £/6 towards — oo, since £ oc |T — T c \~ l l 2 and b is a negative TP constant. 
We distinguish two regions according to the thickness / of the SC sheath near the TP. In 
the first region I = indicating that the stable solution is the null solution without a SC 
sheath near the TP, while in the other region, the stable solution describes a macroscopic 
sheath on both sides of the TP and thus I = oo. Note that a solution with a strictly finite 
thickness, i.e. < / < oo, is never stable and plays no significant role. The two regions 
are separated by the equilibrium surface phase transition line FD (first-order derealization) 
while the dashed line ML marks the metastability limit of the N state of the TP, when T is 
increased towards T c . 

To identify the loci of the first-order phase boundary FD we can apply the technique 
of phase portraits which is a well-known concept in the study of wetting phase transitions 
and has proven to be equally instructive for investigating interfacial phase transitions in 
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superconductors |], |], |16| . In practice, however, it is simpler to calculate 7 from Eq. (^) 
for a macroscopic SC sheath and to see where it equals the surface free energy of the null 
solution. Since 7 = in the latter state the condition for the derealization transition reads 



Itp/sc + Isc/n = 0, (17) 

where 'Jtp/sc is the surface free energy of the bulk SC phase against the TP and ■Jsc/n the 
surface tension of the SC/N interface. The former can be calculated by considering the first 
integral of Eq. ( [TBI) (Jot h = h c ) 

0=7^(1-<A (18) 

which leads, by inserting into the integral (^), using also (0) and considering the half space 
x > 0, to 



2V2 2 
ITP/SC = -T, 



2V2 J 



(19) 



Here, (p + is given by the solution of 

ipl + V2£<p+/b -1=0. (20) 

This equation follows from combining the first integral with the boundary condition (|13|). 
Concerning the surface tension Jsc/n we can use the accurate analytic expression pioneered 
T7\ and improved in ||18|| , 



111 



2 /2 

lsc/N = -X 1.02817^ - 0.13307*;^ + 0(k 2 ^k). (21) 

Due to the high accuracy of this expansion when truncated at order k^Jk in the entire 
type-I regime, we immediately find, by inserting (|l9|) and ([H]) into the condition (|l7|), 
an accurate analytic result for the first-order phase boundary FD. The deviation from the 
numerical results lies within the thickness of the solid line in Fig. |l], even at k = 1/V2. 
For k = the transition occurs at (£/fr)* = —0.6022 and expanding the phase boundary 
about this point reveals that it approaches the k = axis in a parabolic manner k(£/6) ~ 
a {£,/b ~ iC/b)*) 2 , with a ~ 4.95. This demonstrates that in the 1ow-k regime the system 
behaves both qualitatively and quantitatively precisely like the semi-infinite system with a 
wall 1, g |. 

A final aspect of the phase diagram relates to the calculation of the metastability limit ML 
for which it is justified to use the linearized version of the GL theory. In this approximation 



the non-linear terms in the GL equations are omitted so that Eq. (|8|) reads, for T < T c , 

9o = -V?o + al^o/n 2 , (22) 

while the second one becomes trivial with the general solution 

a(x) = ao(x) = h(x + Xq), (23) 

where the boundary condition (|TJ) immediately gives x = 0. Eq. (|23|) expresses that the 
magnetic field is at no point expelled by the SC sheath that nucleates at ML. Thus we need 
only solve Eq. (|22|), subject to the boundary conditions 

0o(O + ) = ^o(0 + ) , <A,(oo) = 0. (24) 

To find (po(x) we follow Ref. J| and reduce fl22| ) to a first-order (non- linear) differential 
equation by introducing the function qo(x) = (p (x) / (p (x) obeying the equation 

go + q 2 = -1 + al/n\ (25) 

with the boundary condition qo(0) = (,/b. This equation must be solved such that qo(x) has 
the acceptable asymptotic behaviour qo(x) ~ —hx/n, implying a Gaussian decay for (po(x). 
This is done by performing (backwardly) the numerical integration of the auxiliary equation 

g + q 2 = -1 + (hx/n) 2 , (26) 

starting from qo(x) = —hx/hi for large x, down to x — 0. For given h/n, the ML value for 
^/b then simply follows from qo{0) = Moreover, using the result for the function qo{x), 
we can construct explicitly a solution to (^) of the form 

<p (x) = <po(0) exp (£ q (x')dx'^ . (27) 

The line ML is obtained by taking a value of k and applying the procedure for the coexistence 
value h c — 1/ \/2. Further, the same method applies away from coexistence (as well as to 
T > T c ) where it serves to obtain the critical nucleation field as a function of temperature 
(see below). 

The phase diagram discussed above demonstrates that for all type-I materials with an 
internal transparent TP there exists an interface derealization transition at some tempera- 
ture Tjj strictly below T c . This transition can be interpreted as a genuine wetting transition 
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such that for temperatures T > Tn the SC state completely wets the TP and this occurs on 
both sides of the TP. For a transparent TP this phase transition is of first order regardless 
of the value of k. This is very different from the analogous phase diagram for wetting at the 
surface of the material for which both first- order and critical transitions are predicted ]TJ, 0. 

A TP displaying TPS is characterized by a negative extrapolation length, hence it is not 
relevant to consider here positive b values. Explicit surface free-energy calculations further 
reveal that for b < the SC phase is preferred by the TP such that ^tp/sc < Itp/n, with 
1tp/n the surface free energy of the N phase in bulk, according to fll4l) . This equality is 
reversed for all k when £/6 > so that the question of wetting by the SC phase is only 
relevant for b < 0. In other words, the line of reversal of preferential adsorption in the 
(k, £/6)-plane is located at £/6 = 0. Also in this respect the wetting phase diagram differs 
from that for an external surface or wall @. 

C. Off-of-coexistence phase behaviour 

Now we turn our attention to the issue of phenomena outside coexistence which can most 
easily be clarified by inspecting the magnetic field versus temperature phase diagram for a 
given material. These diagrams are much more accessible to experimental verification than 
the global diagram in Fig.|l] and provide a means of comparing our results with the TPS 
phase diagrams of Ref . || . 

Fig.||] shows a typical example for k = 0.3 where we employ units based on the dereal- 
ization field Hp = H c (Td) and the derealization temperature Tp, i.e., we take the ratios 
H/ Hd and t = (T — T c ) / (T c — Tp). A first ingredient in this diagram is the bulk coexistence 
line CX with two anchor points, the derealization point D at H/Hp = 1 and t — — 1 
and the bulk critical point at H/Ho = and t = 0. Since the derealization transition is 
first-order the point D is the starting point for a prewetting line FN (first-order nucleation) 
which marks a first-order transition between the null solution with I = and a finite sheath. 
The line FN is tangential to the coexistence line CX at D and changes into a continuous or 
critical nucleation line CN at a tricritical point TCP. The short stretch CN can be computed 
using the technique discussed earlier to determine the line ML in Fig. [l|. It is important to 
stress that this transition is critical since the sheath will appear with an infinitesimal ampli- 
tude. Concerning the spatial extension in the direction perpendicular to the TP, however, 
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it turns out that the sheath always has a thickness of the order of £ even just below the 
transition. Note that at the tricritical point FN and CN meet with common tangents and, 
to the left of TCP, CN continues as the metastability limit ML of the N state of the TP 
(dashed line). Finally, the line CN ends at zero field at the temperature T c tp which obeys 
the equation |TJ, |2j 

cc^tp) = _ 1; (28) 

showing that the extrapolation length b can be found in principle from experimental deter- 
mination of T C] tp- 

We obtain the first-order nucleation line FN by numerically solving the GL equations for 
a sheath-type solution and finding the point in the phase diagram at which the free energy 
of this solution is zero. The location of TCP can be determined accurately by extending 
the zeroth-order solution of the linear GL equations (see above) a little further. In order to 
obtain a non- vanishing solution of the non-linear theory we need a correction to ao(x) and 
<Po(x). In view of Eq. (|27f) <po(x) is the small quantity through its small amplitude and thus 

it suffices to expand a(x) = a (x) +ax(x) + Using the first integral of the equations and 

working to second order in ip we immediately get 

ai = ^(^-*8±<P§), (29) 



or, using 



o-i = ^qo(x)<fl(x). (30) 



This result is particularly useful when we use it in an alternative expression for 7 evaluated 
in the extrema (see 0) 

7 = /_^dz{-|- + (d-/i) 2 }. (31) 

The advantage of this expression is twofold, firstly the boundary term is absorbed and 
secondly it shows that both terms are small, of order <£>q for small (po. Along the line CN, 
7 = because <fo = and do = h. Another mechanism, however, to make 7 vanish is the 
compensation of the two terms in (|3"T|), which can be put in the form 



/+00 r+00 
dxq 2 (x)v 4 (x)/ / dxipt(x). (32) 
-00 J —00 

and which is of use in determining TCP as follows. For each value of k a continuous set 
of (h,£/b) pairs can be found from the linear theory, but only one pair will satisfy the 
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condition (|32]) for the onset of the nucleation of a finite (in contrast with infinitesimal) 
sheath, which is then by definition the TCP. 

The nucleation transitions discussed above indicate the point at which local supercon- 
ductivity sets in near the TP, which obviously corresponds to the phenomenon of TPS. 
Consequently the results presented here should agree with the predicted TPS phase dia- 
grams for a parallel magnetic field, at least with those pertaining to strongly transparent 
TP's f|. We have checked this for various k values and indeed found an excellent agree- 
ment for all the transitions. Furthermore, we may interpret the TPS transition as a genuine 
prewetting transition and thus by reinterpreting the existing theoretical and experimental 
TPS phase diagrams we stress that prewetting has since long been observed in classical 
superconductors. Intriguingly, we recall that only in more recent years has clear evidence 



of prewetting transitions been found in experiments on classical binary liquid mixtures [19 



An important feature of the prewetting line found here is that it does not end at a surface 



critical point as typically predicted and found in liquid mixtures [|1J, [L9| but becomes, via a 
tricritical point, a continuous nucleation line. Further, our results are particularly interest- 
ing with respect to the behaviour near the wetting point D, since from the experimental TPS 
diagram it is unclear what happens to the TPS transition near the bulk critical field H c . 
We now understand that the surface phase transition at the TP becomes a bulk transition 
precisely at the derealization transition. 



IV. OPAQUE TWINNING PLANES 

A. Classification of solutions and terminology 

In this section we focus on the second limiting case, i.e., that of completely opaque TP's 
for which axp — > oo. In this case the two sides of the plane are, to some extent, independent 
since the requirement of a continuous order parameter at x = is no longer applicable, 
hence for a general solution (p + ^ ip_. However, there remains some coupling between the 
two half-spaces due to the vector potential and the magnetic induction which must both be 
continuous at the TP. Moreover, the number of possible solutions is further restricted by the 
assumption that the extrapolation length b is equal on both sides of the TP. This means that 
the critical temperature T Cj tp (which is defined for zero magnetic field) is assumed to be the 
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same for both sides. Nevertheless, a large number of different solutions can still be found 
among which several are characteristic for the opaque limit and deserve special attention 
here. For a neatly arranged description of the phases and phase transitions it is appropriate 
at this stage to present a classification and a schematic overview of the various solutions at 
bulk coexistence with equal bulk (N) conditions on the two sides of the TP. 

A first (trivial) solution is the null solution <p(x) = without local superconductivity. 
For solutions with at least one finite SC sheath we distinguish three scenarios as depicted 
in Fig. |3|. Interestingly, due to the presumed opacity we can consider imposing <p(x) = for 
one half-space, say x < 0, leading to a solution as shown in Fig. |](a) and referred to as "wall 
solution" since it corresponds exactly to the solution found in a system with a wall obeying 
the boundary condition d(0) = h. In the case of a double SC sheath we can have either 
a continuous solution with tp + = and a(0) = (Fig. 0(b)) corresponding to the sheath 
solution for a transparent TP, or a more general discontinuous solution with ip + ^ ip_ and 
a and a continuous at x = (Fig. 0(c)). Arguing along the same lines we also find a series 
of solutions with a macroscopic SC layer on only one or on both sides of the TP. For the 
former, in the other half-space, we can either have the wall solution with ip(x) = or a finite 
sheath with a continuous or discontinuous order parameter at the TP. The order parameter 
profiles for these solutions are schematically drawn in Fig. [|. For a double macroscopic 
layer, on the other hand, we should in principle take into account both the continuous and 
discontinuous solution. It can be shown, however, by a simple argument (see below) that 
the latter is impossible and therefore unnecessary to consider. As outlined in the previous 
section a double macroscopic layer will completely expel the magnetic field such that a = 
over the entire region of the SC layer. Consequently expression (|l8f) applies also here and 
at x — +, ~ this yields 

0+,- = ^(l-<_), (33) 
While it may appear that any combination of ip + and ip_ obeying these equations is a 
possible solution, this is only true if the boundary conditions at the TP are omitted. Indeed, 
by applying the boundary conditions which in this case are again given by Eq. (|13|) it is 
straightforward to see that, since b is assumed to take the same value on both sides of the 
TP, ip + must equal for a solution with two macroscopic SC layers as exemplified in Fig. |5|. 
Note that this solution is the one we found for the case of a transparent TP. 

From the above classification of solutions it is obvious that an opaque TP is in a sense a 
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combination of a system with a wall and one with a transparent TP, a feature which gives 
rise to additional phase transitions. Before further embarking on the phase behaviour it is 
convenient to comment here briefly on the terminology used below. The transition from the 
null solution to a state with either one or two finite sheaths as well as the transition from 
a single sheath to a pair of sheaths is referred to as a nucleation phenomenon. Likewise, 
the transition from a "wall state" with one macroscopic SC layer to a state that is the 
combination of such a layer and a finite sheath is called nucleation. Further, the appearance 
of one macroscopic SC layer, either from the null solution or from one finite sheath, is 
a delocalization or wetting transition, while a depinning transition refers to the transition 
between one and two macroscopic layers. The latter describes the depinning of a SC/N 
interface that is initially pinned at the TP and is analogous to the one investigated in 
Ref. [Q. We start with a discussion of the TP states at bulk two-phase coexistence. 

B. Phase diagram at bulk coexistence 

The fundamental diagram for interfacial phase transitions in the case of opaque TP's is 
presented in Fig. 0. The different regions are distinguished according to the values of the 
two length scales l\ and I2 representing the thicknesses of the sheaths on the two sides of the 
TP. We assume here, without loss of generality, that l\ > 12- It is striking how different this 
phase diagram is compared to that for the case of perfect transparency, see Fig. p]. As we 
will demonstrate below, the main part of the diagram for the opaque case simply recovers 
the results of the system with a wall It is important to stress at this stage that we are 
only interested in the wall results for negative b values since this is the relevant region for 
TP's that show TPS. 

One ingredient of Fig. concerns the stability of the sheath solutions given in Fig. |3|. 
From our analysis in Section [TTT] we know that at bulk coexistence the double symmetric 
sheath (Fig. 0(b)) is never stable, and hence plays no significant role. Furthermore, explicit 
free energy calculations reveal that the discontinuous double sheath (Fig. 0(c)) has an even 
higher free energy than the continuous one and thus is irrelevant. The single sheath solution 
(Fig. |3|( a)) is stable for a certain interval of £/& values provided that k > 0.374 as can be 
inferred from the results for a system with a wall [p], g|. Hence, in the upper right corner of 
the diagram we retrieve a feature of the wall diagram, with a finite sheath on only one side 
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of the TP while the other side is still in the N phase, thus h > and I2 = 0. In the lower 
half of the diagram we have on the low-£/|6| side a region where no stable sheath solution 
exists other than the null solution (i x = and I2 = 0). The two regions are separated either 
by the line CN (critical nucleation) ending in a tricritical point TCP, or by the short stretch 
FN (first-order nucleation) between TCP and a critical end point CEP. To the left of TCP 
the line CN continues as the metastability limit MLi of the null solution. 

Next we need to consider the solutions of Fig. |] and f| and to investigate their stabil- 
ity. Our calculations suggest that when the temperature is increased towards T c the system 
first enters the state given in Fig. f§(a), either from one finite sheath or from the null so- 
lution. These transitions are indicated by the line CD (critical derealization) and FD 
(first-order derealization) respectively, meeting each other at the CEP. By further increas- 
ing the temperature the system undergoes a first-order depinning transition towards the 
double symmetric solution, i.e. a transition from the solution given in Fig. |](a) towards the 
solution of Fig. |5[ Thus the regions li — 00, I2 — and l\ = I2 — 00 are separated by the 
first-order phase boundary FDP (first-order depinning). Lastly, the line ML2 represents the 
metastability limit of the state with 1% = 00, I2 = 0, with the SC/N interface pinned at the 
TP. Note that the solutions represented in Fig. |](b) and (c) are never stable. 

An important conclusion that we draw from this phase diagram is that for relatively 
low temperatures an opaque TP behaves as a system with a wall, with on one side of the 
TP the different wall solutions while the other side remains fully in the N phase. In other 
words, to the right of the line FDP in Fig. ^ we immediately obtain the phase diagram by 
copying the results for a wall without additional computations. All the details concerning 
the determination of the transition lines for this part of the diagram can be found in Ref. @ 
and ||. In particular, the critical derealization condition reads 

Iw/n = Iw/sc + Isc/n, (34) 
while the first-order derealization condition is simply given by 

Iw/sc + Isc/n = 0, (35) 

where jw/n (jw/sc) is the surface free energy of the bulk N (SC) phase against a wall. The 
novel feature in Fig. || is the FDP line arising from the possibility of having a macroscopic 
SC layer also in the other region leading to a depinning transition at higher temperatures 
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along with the associated spinodal ML2. It is precisely these two transition lines that we 
will be concentrating on in the remainder of this section. 

The condition for the first-order depinning transition is obtained by equating the free 
energy of the solution with one macroscopic layer at x > and no sheath at x < to that 
of the solution with the double symmetric macroscopic SC layer. The former can be written 
conveniently as the sum jw/sc + lsc/Ni thus the depinning condition reads 

Iw/sc + 1SC/N = ^{lTP/SC + 1sc/n)- (36) 

In general the depinning phase boundary has to be calculated numerically. An approximate 
analytic result can be obtained using the powerful expansions in k for the surface free 
energies. For the surface tension ^sc/n we use the result ( pT|) while a similar expansion has 
been derived for jw/sc m Ref- 0- Using these earlier results in the depinning condition (|36|) 
provides a very accurate approximation for the phase boundary across the complete type-I 
range with an error less than 1%. For k = the depinning and the derealization transitions 
coincide at (£/&)* = —0.6022. Both phase boundaries have a parabolic foot near k — 0, i.e. 

~ a (£/fr — (£/fr)*) 2 with a 27.0 for depinning and a ps 4.95 for derealization |5J. 
From a mathematical point of view it is interesting to examine whether the extensions of 
the lines FDP and CD meet in the type-II regime, for k > 1/ \[2 . The condition for a line 
crossing is given by combining fl3*3|) and (^). Explicit calculations using an expansion in k 
for the line CD obtained in reveal indeed an intersection at k ~ 0.815, £/6 ~ —0.251. 

The line ML 2 marks the nucleation of an infinitesimal sheath on one side of the TP under 
the condition that on the other side of the TP a macroscopic sheath is stable. To obtain 
this line we first compute numerically a solution for a macroscopic layer in one half space, 
x > say, for a given value of £/& subject to the condition d(0) = h. To proceed we combine 
it with a solution of the linear GL theory applied to the region x < using the scheme 



introduced in Section |T|. The zeroth-order solution for the vector potential is again given 
by Eq. (|23| ) with in this case Xq determined by the value of a(0) which we get from the 
numerical calculation in x > 0. Finally we have to compare ^0(^0) with the given value 
of £/& and iterate the procedure such that the two become equal. This defines the critical 
nucleation of the sheath in the region x < 0. The line ML2 is now formed by applying 
this method for the coexistence field h c . The procedure is equally applicable for any other 
field appropriate for the study of off-of-coexistence phenomena, an issue which is addressed 
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below. 



C. TPS phase diagrams 

From the properties of an opaque TP at bulk coexistence elucidated above we anticipate 
that the onset of local superconductivity near this plane generally happens in two steps. This 
is a natural consequence of the stability of "wall solutions" in these systems with <p(x) = 
on one side of the TP while on the other side the SC phase has already nucleated. We now 
wish to concentrate on the various nucleation phenomena, the order of which depends on 
the parameter Kj^ clS 8b function of magnetic field and temperature. 

For k < 0.374 the derealization transition on one side of the TP is first-order and thus 
will have a prewetting extension into the region of the phase diagram where the N phase 
is stable in bulk. This prewetting line FN is the first-order nucleation of a sheath on one 
side of the TP, and hence corresponds to the transition between the null solution and a 
solution as given in Fig. 0(a). This line will be tangential to the bulk coexistence line 
CX at the first-order derealization point D and changes at the tricritical point TCP into a 
continuous nucleation line CN. The latter describes the nucleation of one infinitesimal sheath 
and continues as the metastability limit MLx of the null solution to the left of TCP. So far 
the surface phase diagram, an example of which is given in Fig. [7] for k = 0.3, resembles the 
diagram for the system with a wall [l], 0] . Note that we use the same units as in Fig. 0. 

As we have argued above, however, another distinct point exists at bulk coexistence at 
higher temperatures but still below T c , namely the first-order depinning point DP. The 
presence of the depinning transition leads to a variety of additional transitions compared 
to the wall system. Indeed, due to the first-order character of this point a second first- 
order nucleation line FN 2 appears, attached to the point DP. This phase boundary can be 
interpreted as a predepinning line. For the sake of clarity this line is omitted in the main 
figure, with all the details near the depinning transition given in the inset of Fig. [T[ At 
the line FN 2 the SC phase nucleates on the other side of the TP, hence it is the transition 
between one sheath and a double sheath. The line FN 2 is tangential to the line CX at DP 
and meets a continuous nucleation line CN 2 at a second tricritical point TCP'. To the left of 
TCP', CN 2 continues as the metastability limit ML 2 of the state with one sheath. Obviously 
CN 2 denotes the critical nucleation of the second SC sheath. The two critical nucleation 
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lines end in zero field at the same point i C) t p which follows from the assumption that the TP 
is characterized by a single value of the parameter b. 

For materials with k > 0.374 for which the derealization transition is critical, the field- 
temperature phase diagram will undergo qualitative changes as regards the wetting-like 
transitions, i.e., the transitions related to the nucleation of the SC phase on one side of the 
TP. We refer interested readers to Rcf. || for examples of (H, T)-phase diagrams containing 
these nucleation transitions and restrict ourselves here to a brief discussion of the pertinent 
features relevant for the present study. For k = 0.374 the point D and the starting point 
of the first-order nucleation line FN separate and by further increasing k this first-order 
transition disappears. The critical nucleation line CN then extends to low temperatures 
without intersection with the bulk coexistence line CX. Hence in this case the nucleation 
of the first sheath is always critical. For the nucleation of the second sheath there is still 
the possibility of first-order as well as second-order transitions. In fact, we know that 
the depinning transition is first-order for all type-I materials so that the situation that we 
sketched for k = 0.3 near the depinning point DP is representative for the entire type-I 
regime. Thus for an opaque TP there are always two distinct nucleation transitions, i.e. the 
SC phase never appears simultaneously on both sides of the TP. 

V. DISCUSSION AND CONCLUDING REMARKS 

In this paper we have analyzed the phase behaviour of the SC/N interface near internal 
TP's in type-I superconductors. Our calculations reveal that the results are highly sensitive 
to the boundary conditions imposed at the TP, which in turn depend on the degree of 
transparency of the TP for electrons. For perfectly transparent TP's the order parameter 
is continuous such that only symmetric profiles need to be considered and at first sight 
the analysis resembles that of an external surface or wall. The only differences originate 
from the boundary condition for the vector potential and, surprisingly, this small technical 
modification turns out to have dramatic consequences for the order of the derealization 
transition. In particular, for a transparent TP only first-order transitions are found for the 
entire type-I regime while the possibility of critical transitions that are predicted for the 
wall system |], |2j are suppressed. Our results for the magnetic field versus temperature 
phase diagrams are in excellent agreement with earlier experimental and theoretical results 
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obtained in the context of TPS ||. 

We comment further that our approach for the nucleation of the SC phase near the TP 
equally applies to type-II materials (with k > l/y/2). Specifically, we have considered the 
case of Nb with k ~ 1 for which the experimental results suggest that the properties of the 
TP are very close to the limiting case of perfect transparency. The nucleation transition in 
this case is always of second-order and our results again perfectly agree with earlier work || . 
We remark that in this case the nucleation lines are only relevant at sufficiently high fields 
above their intersection with the upper critical field H c2 . 

The situation drastically changes when considering opaque TP's in which case a discon- 
tinuity in the order parameter is allowed at the TP. The decoupling of the two sides of 
the TP makes it possible to consider wall solutions with tp(x) = on one side of the TP. 
Moreover, from free-energy considerations, we have shown that these solutions are stable 
in a large region of the phase space and, as a result, the system precisely undergoes the 
transitions predicted for a wall system. In this case only one side of the TP will be wetted 
by the SC phase at the derealization transition which can be either first-order or critical 
depending on the value of k. By further increasing the temperature at bulk coexistence a 
first-order depinning transition is predicted for all type-I materials. Consequently, the field- 
temperature phase diagrams for opaque TP's fundamentally differ from their counterparts 
for the transparent limit. A characteristic feature of an opaque TP is the existence of two 
distinct nucleation lines which in principle should be measurable and thus can provide an 
experimental means of distinguishing between transparent and opaque TP's. Related to this 
we remark that experiments in Sn (k ~ 0.13) have revealed only one nucleation transition, 
although it is assumed that the twin boundary in this material has a low transmission for 
electrons [§, 12], [U| . The apparent absence of a second nucleation transition can in this case 
be attributed to the low-K value of the material, since in the 1ow-k regime the various tran- 
sition lines lie extremely close together (see Fig. H) and they would be difficult to distinguish 
in an experiment. Note that in the limit k —>■ the differences between a transparent TP, 
an opaque TP and an external surface vanish. 

Finally, for type-II materials with opaque TP's our calculations show that the nucleation 
can be either first-order or critical and this is at variance with the conclusion of earlier 
work [BT| stating that the TPS transition is always second-order for type-II materials. Our 
study demonstrates that this surmise is correct only for the case of transparent TP's. The 
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reason for this is that the tricritical nucleation point TCP merges with the derealization 
transition at k — l/v2 for transparent TP's. This can be seen also from the merging of the 
spinodal ML with the line FD in Fig. [I]. In contrast, for opaque TP's the TCP of nucleation 
remains well off of coexistence, at H > H c , and in the type-II regime at H > H C 2, so that 
there is room for first-order nucleation. This is further demonstrated by the fact that, in 
Fig. [| the lines ML 2 and FD are still far apart at k = 1/ \/2. 

One remarkable implication of our work is that, in general, with the exception of perfectly 
transparent TP's, there exist stable states of local superconductivity which are asymmetric 
about the TP. The possibility that a SC sheath or a macroscopic SC layer can occur on 
one side of the TP while the other side is in the normal state is indeed noteworthy and 
has been met with scepticism. It has been suggested that, since our analysis is essentially 
one-dimensional and neglects states which are inhomogeneous in the direction(s) parallel to 
the TP, there may exist modulated states , e.g., composed of a linear array of soft vortices 



parallel to the TP [ 1 3 1 or states with a local field penetration and a change of phase of the 



wave function, which may have a lower free energy than the states we have considered [22 



Although we cannot rule out this possibility at present, we would like to point out that 
asymmetric wetting, followed by symmetric depinning, has been found previously in the 
context of grain-boundary wetting in the framework of a real scalar order parameter 
model of Ising type. 
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FIG. 1: Wetting phase diagram at bulk two-phase coexistence for a transparent twinning plane 
(TP) in the variables k (with < k < l/\/2) and £/b (with b < 0). The solid line represents the 
first-order derealization transition and separates the region without a SC sheath (I = 0) from that 
with a macroscopically thick SC layer (I = oo). The dashed line is the metastability limit of the 
null solution corresponding to the normal state of the TP, encountered when, e.g., the temperature 
is increased towards T c . 
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FIG. 2: Magnetic field versus temperature phase diagram for a transparent TP with k = 0.3 in 
the variables H/Hn and t = (T - T C )/(T C — To). At the nucleation transitions FN and CN local 
superconductivity appears at the TP. Further details of the phase boundaries are provided in the 
text. 
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FIG. 3: Sketch of the order parameter profile for sheath- type solutions: (a) Wall solution with 
<p(x) = for x < and d(0) = h. (b) Continuous solution with ip + = (p_ and a(0) = 0. (c) 
Discontinuous solution with ip + / ip- and a and a continuous at x = 0. 
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FIG. 4: Sketch of the order parameter profile for solutions with one macroscopically thick SC layer: 
(a) Wall solution with <p(x) = for x < 0. (b) Continuous solution with ip + = and a finite 
sheath for x < 0.(c) Discontinuous solution with ip+ / ip- and as in (b) there is a finite sheath for 
x < 0. 25 
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FIG. 5: Schematic order parameter profile for a continuous solution with a double macroscopic 
superconducting layer. Both sides of the TP are in this case wetted by the SC phase. 
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FIG. 6: Wetting phase diagram at bulk two-phase coexistence for an opaque TP in the variables 
k and £/b. The different phases are characterized by the values l\ and I2 of the thicknesses of the 
SC sheaths on the two sides of the TP. The various transitions and transition lines are explained 
in the main text. 
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FIG. 7: Surface phase diagram for an opaque TP with n = 0.3 in the variables H/Hn and 
t = (T — T C )/(T C — Tr>). The inset shows the different nucleation transitions near the first-order 
depinning point DP. All the transitions are explained in the text. 
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